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Abstract 

At NNLO it is particularly important to have a Variable-Flavour Number Scheme 
(VFNS) to deal with heavy quarks because there are major problems with both the zero 
mass variable-flavour number scheme and the fixed-flavour number scheme. I illustrate 
these problems and present a general formulation of a Variable-Flavour Number Scheme 
(VFNS) for heavy quarks that is explicitly implemented up to NNLO in the strong cou- 
pling constant as, and may be used in NNLO global fits for parton distributions. The 
procedure combines elements of the ACOT(x) scheme and the Thorne- Roberts scheme. 
Despite the fact that at NNLO the parton distributions are discontinuous as one changes 
the number of active quark flavours, all physical quantities are continuous at flavour 
transitions and the comparison with data is successful. 



1 Introduction 



While up, down and strange quarks may be treated as being effectively massless partons, the 
heavy quarks, charm ~ 1.5GeV, bottom ~ 4.3GeV, top ~ 175GeV, must have their mass, 
tuh taken into account in any QCD calculations. In particular it is essential to treat charm 
and bottom correctly in global fits for parton distributions. There are two distinct regimes 
that can be considered. Near threshold, Q 2 ~ m^, massive quarks are not treated as parton 
constituents of the proton but are created in the final state. Any processes may be described 
using the Fixed-Flavour Number Scheme (FFNS). For example, structure functions are given 
by 

F(x,Q 2 ) = C riQ'/m^Q 2 /^)®^^), (1) 
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up to higher twist (0(A 2 /Q 2 )) corrections, where rif is the number of light partons and all the 
mass dependence is in the hard coefficient functions which have been calculated up to NLO 
(i.e. 0(a 2 s )) [1]. This is reliable for scales not much greater than m 2 H , but increasing orders 
in as contain increasing logarithms in Q 2 /m 2 H , and order-by- order perturbation theory is not 
guaranteed to be accurate. Also, the FFNS coefficient functions are not known yet at NNLO, 
rendering an NNLO FFNS impossible to define. 

At high scales, i.e. Q 2 ,fi 2 ^> rrijj, the heavy quarks are expected to behave like massless 
partons. The heavy quark is treated like the other partons and ln(/z 2 /m^) terms are then 
automatically summed via evolution. The simplest approach is the Zero Mass Variable-Flavour 
Number Scheme (ZMVFNS) [2]. This ignores all 0(m 2 H /Q 2 ) corrections for each of the rin 
heavy quarks , and the structure functions are given by 

F(x,Q 2 ) = C? MVF (Q 2 /» 2 ) ® /;' + "'V), (2) 

where the hard coefficient functions are mass independent. Although this is called a "scheme" 
it is important to note that unlike usual scheme definitions, which are alternative ways to order 
the perturbative series, the ZMVFNS is incorrect by terms of 0(m 2 H /Q 2 ), and is really only an 
approximation in the region m 2 H ~ Q 2 . The approximation in this region may indeed be very 
important in practice, as I will demonstrate later. A correct variable-flavour number scheme 
should not have these inaccuracies, but should correct the coefficient functions for the mass 
effects. 

As we go from an n/-flavour to an rif + 1-flavour scheme, the partons in the different number 
regions are related to each other perturbatively, 

f: i+ \f) = A 3k ^ 2 /m 2 H )®f;i^% (3) 

where the perturbative matrix elements Ajk{[i 2 / m 2 H ) contain the \n(ii 2 / m 2 H ) terms which relate 
fk f {^ 2 ) an d /fc /+1 (A* 2 ) an d lead to the correct evolution for both. There is then a similar 
relationship as we go from ann/ + 1-flavour to an rif + 2-flavour scheme, but I will consider the 
transitions one at a time in this paper. At the bottom quark transition point rif is effectively 
equal to 4, i.e. the charm quark is already evolving like a massless partons below this point. 

At LO, i.e. zeroth order in as, the relationship is trivial, 

g^rV) = ^):v). (4) 

At NLO, i.e. first order in as, the non-trivial contributions are 

(h + h)(v 2 ) = ^P° g ®g n f(Q 2 )\n(v 2 /m 2 H ), ff "' + V) = (l - ^ ln^/m 2 ^ V), (5) 

where h(x, fi 2 ) is the heavy quark parton distribution. Hence, the heavy flavour evolves from 
zero at fi 2 = m 2 H according to standard massless quark evolution and the gluon loses correspond- 
ing momentum. It is natural to choose fi 2 = m 2 H as the transition point from the n/-flavour to 
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Figure 1: The discontinuity in F 2 c (a;,Q 2 ) (left) and F 2 (x,Q 2 ) (right) using the zero-mass 
variable- flavour number scheme at NNLO. 

the rif + 1 -flavour scheme, since at this order the partons are then continuous and the heavy 
quark starts evolving from a zero value. At NNLO, i.e. second order in as, there is much more 
complication 

f^NLo^ 2 ) = (g) 2 E(4° + 4 1 in(/i 2 /<) + 4 2 ^ 2 (^/<))®/; / (/i 2 ) ) (6) 

ij 

where the A 2 j° [3] are generally nonzero. There is no longer any possibility of a smooth transition 
at fi 2 = m 2 H . Since each of the A^ is a different function there is no smooth turn on of the 
flavour distribution at any alternative value of fi 2 , hence, for technical simplicity, it seems 
sensible to keep the transition point at /i 2 = m 2 H , though any other point can be chosen 
if desired. Making this choice it turns out that A 2 jf g , the matrix element giving the gluon 
contribution to the heavy quark distribution, is negative at small x, even though the structure 
function is always positive, and the heavy quark starts evolving from a negative value in the MS 
scheme. This highlights the fact that parton distributions are not physical quantities. However, 
it also illustrates a major problem with the zero-mass variable-flavour number scheme. 

In order to make a concrete illustration of the effect we must choose a factorization and 
renormalization scale. For light partons both of these are conventionally chosen to be /i 2 = Q 2 . 
It is most natural to place the heavy flavour on the same footing, and choose the same scale, 
e.g. it is difficult to think of momentum conservation at a given scale if the factorization scale 
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is different for different partons. Hence, for the remainder of this article I will set /i 2 = Q 2 . 
Having made a choice, we see that in the ZMVFNS the coefficient functions already lead to a 
discontinuity in the structure functions at NLO, i.e. 

F 2 H (x,Q 2 ) = Q 2 <m 2 H , 

= ^C 2 , g ®g n f +1 (Q 2 ) Q 2 >m 2 H . (7) 

However, this is a very small effect at NLO. It is larger at NNLO, since the coefficient function 
makes a larger contribution, even though it is nominally of higher order, and it is negative at 
smallish x (x ~ 0.001), as is the starting value of the partons. Hence, F 2 H (x, Q 2 ) is also negative 
for Q 2 = m 2 H . This shows that the ZMVFNS is not really feasible at NNLO. It leads to huge 
discontinuities in the charm structure function F 2 (x, Q 2 ), and even significant discontinuities 
in the total structure function F 2 (x, Q 2 ). These are shown in Fig.l, where a small discontinuity 
in F 2 (x, Q 2 ) at Q 2 = m 2 is also seen. This is a measure of the mistake made in omitting 
the 0(m 2 H /Q 2 ) corrections in this approximate "scheme". One really needs a general Variable- 
Flavour Number Scheme (VFNS) joining the two well-defined limits of Q 2 < m 2 H and Q 2 3> m 2 H 
in a theoretically correct manner. We will outline such a scheme, which has been implemented 
explicitly up to NNLO, in the remainder of this article. 2 



2 The Variable-Flavour Number Scheme 

A correct VFNS can be defined by demanding equivalence of the rif (FFNS) and rif + 1-flavour 
descriptions at all orders, i.e. 

F(x,Q 2 ) = Cr(Q 2 /m 2 H ,Q 2 /^ 2 ) ® #V) = Cj F (Q 2 /m 2 H ,Q 2 /^ 2 ) ® /?' + V) 

ee CY{Q 2 lm 2 H M 2 l^)®A^ 2 lm 2 H )®f k \ l i 2 ). (8) 

Hence, the VFNS coefficient functions have to satisfy 

CF(Cf/m 2 H, OV) = Cj F (Q 2 /m 2 H , Q 2 /^) ® A jk (^/m 2 H ), (9) 

at all orders. It is important to remember that the left-hand side of this expression is expanded 
in the n/-flavour coupling ag jTV/ , while the right-hand side is most naturally expanded in the 
rif + 1-flavour coupling cxs,n f +i- The two are related by 

«5,n /+ l(/i 2 ) = a S ,n f + H» 2 /™%) <X%n f + 0(<*S,n f )- ( 10 ) 

The coupling is therefore continuous up to 0{a 2 s ), but at 0{a z s ) there is a small discontinuity. 
This discontinuity does not influence the VFNS up to NNLO. 

2 The scheme has previously been outlined in a very brief form in [4]. 
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At 0(a s ) eq.(9) becomes, for example, for the structure function F 2 (x,Q 2 ) 

C^\Q 2 /m 2 H ) = Cl F / H {Q-lm\)®Pl g \^lm^ (11) 

The VFNS coefficient functions automatically tend to the massless limits as Q 2 /m 2 H — > oo [5] 
and, if we use the zeroth order cross-section for photon-heavy quark scattering, 

CZh'h(Q 2 /<, *) = (! + *m 2 H /Qr)6(z - 1/(1 + m 2 H /Q 2 )), (12) 

this is the original ACOT scheme [6]. 

However, C^nniQ 2 1 ' m ?r) ( we se ^ f^ 2 — Q 2 explicitly for simplicity) is only uniquely defined 
in the massless limit Q 2 /m 2 H — > oo. One can swap 0(m 2 H /Q 2 ) terms between C^'° H (Q 2 /m 2 H ) 
and C 2 ^ 1 [Q 2 / m 2 H ) while still satisfying Eq.(ll), i.e. C 2 ^(Q 2 /ra 2 H ) is not uniquely defined. 
This is true for all C^nniQ 2 1 ' m #) • The original ACOT prescription removed the ambiguity by 
defining C 2 ^' l H (Q 2 /rrijj) as the calculated coefficient function for an incoming massive quark. 
However, this violates the physical production threshold W 2 > 4m 2 H since it only needs one 
quark in the final state rather than a quark- ant iquark pair. Hence, there is not a smooth 
transition at Q 2 = m 2 H asw/ - ► n/ + 1. The Thorne- Roberts Variable-Flavour Number Scheme 
(TR-VFNS) first recognized this ambiguity in the definition of C 2 hh(Q 2 7 m if) 17] an d removed 
it by imposition of the physically motivated constraint of (dF 2 /d\iaQ 2 ) being continuous at the 
transition point Q 2 = m 2 H (in the gluon sector). Hence, it guaranteed smoothness at Q 2 = m 2 H , 
but the approach to Q 2 /m 2 H — > oo is a little odd - the VFNS result overshooting the zero mass 
result before approaching it asymptotically from above. This effect diminishes at higher orders 
but more of a problem is the complicated form of the scheme - C%h%{Q 2 /m 2 H ) oc (P° ) -1 , 
which is not a simple function. This makes the scheme very involved at higher orders and it is 
also not well suited to charged currents [8]. 

There have been various other alternatives since this. Most recently Tung, Kretzer and 
Schmidt have devised the ACOT(x) prescription [9] which may be interpreted as 

C^° H (Q 2 /m 2 H ,z) = 5(z-Q 2 /(Q 2 + 4m 2 H )), 
-> F 2 H '°(x,Q 2 ) = (h + h)(x/x rnax ,Q 2 ), x max = Q 2 /(Q 2 + 4m 2 H ). (13) 

Hence, the zeroth-order coefficient function tends to the standard C 2 hh{ z ) = ^(1 — z ) f° r 
Q 2 /m 2 H — > oo but respects the threshold requirement W 2 = Q 2 (l — x)/x > 4m 2 H for quark- 
antiquark production. Moreover, it is very simple. For the VFNS to remain simple (and 
physically motivated) at all orders n in as it is necessary to choose 

C 2 ,HH(Q 2 / m Hi z ) — C 2 ^u(z/ x max ). (14) 
It is also important to choose 

Cljih{Q 2 I^Hi z ) Cl%h{ z I x max) i (15) 
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and to impose the condition that C^^f H (Q 2 /m 2 H , z) = (as is done in [7]), despite the fact that 
one obtains 

2 

Tn 

Cl HH (Q 2 /m 2 H , z) = 4z^6(z - 1/(1 + m 2 H /Q 2 )) (16) 

for single quark-photon scattering. Because the heavy- flavour contribution to F L (x,Q 2 ) is 
highly suppressed at low values of Q 2 , compared to that for F 2 (x,Q 2 ), we choose 

C I!hh(Q 2 ''/ m H> Z ) = 4(1 _|_4 m 2^/Q2 ~~ -^j C L%H{ Z / X raax)- (17) 

The prefactor is zero at Q 2 = m 2 H and tends to unity as Q 2 /m 2 H — * oo, but is suppressed by the 
physical threshold of 4m 2 H for intermediate values of Q 2 . This factor guarantees that the heavy- 
flavour parton contribution to F^(x,Q 2 ) is heavily moderated for Q 2 just above m 2 H . Failure 
to do this results in a kink in F^(x, Q 2 ) just above Q 2 = m 2 H . The fact that this is a function 
of Q 2 /m 2 H only avoids any additional complications when obtaining higher-order coefficient 
functions by convolutions with matrix elements which the definition in Eq.(9) requires. 

Adopting this convention for the heavy-flavour coefficient functions at NNLO we have, for 
example, 

C2,Hg(Q 2 / m H, z) = C^g(Q 2 /m 2 H , z) - C^hh (z/x max ) <g> A^Q 2 'm 2 H ) 
-C™£(z/x max ) ® A 2 Hg {Q 2 /m 2 H ) - i- \n(Q 2 /m 2 H ) C™£(z/x max ) ® A l Hg {Q 2 /m 2 H ). (18) 

The last term comes from the change in the coupling constant as we go across the transition 
point, i.e. from Eq.(lO). This would be absent if we used (somewhat unnaturally) the n/-flavour 
renormalization scheme above Q 2 = m 2 H) as is sometimes done, but this means the definition of 
A 2 Hg {Q 2 /m 2 H) z) is different in the two renormalization schemes (compare that in [3] with that 
in [10]). There is also in principle a contribution of the form 

~ MQ 2 /m 2 H ) Cl F / g {Q 2 /m 2 H , z) - C^ g (Q 2 /m 2 H , z) ® A] g {Q 2 /m 2 H ) (19) 

on the right-hand side, but 

A l gg (Q 2 /m 2 H ,z) = -±-\n(Q 2 /m 2 H )5(l-z), (20) 

as seen in Eq.(5), so these terms cancel. Both would be absent if we used the n/-flavour 
renormalization scheme above Q 2 = m 2 H . 

From the definition in Eq.(18) we see that since A 2 Hg (l,z) ^ 0, the coefficient function 
C 2 Hg (Q 2 /m 2 H) z) is discontinuous as we go across Q 2 = m 2 H . This compensates exactly for 
the 0{a 2 s ) discontinuity arising from that in the heavy-flavour parton distribution, i.e. for the 
term C^hh ®(h + h), and F^(x, Q 2 ) is continuous. 3 In practice this requires the knowledge of 

3 In principle there are 0(a%) discontinuities due to terms such as C^hh ® {h + h) and C^'g ® g nf+1 , 
i.e. 0(as) coefficient functions convoluted with 0{or s ) discontinuities in partons. These would be cancelled at 
NNNLO by discontinuities in O(o.g) coefficient functions, but are actually tiny in practice. 
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C F jj g (Q 2 /m^, z )- An expression for this exists as semi-analytic code [11] where the dominant 
contributions for W 2 — > oo and W 2 — > 4m 2 H (from above) are analytic and the rest numerical. 
I have produced much faster analytic expressions which are exact for Q 2 /m 2 H — > oo and in 
some cases for W 2 — > Am 2 H) and the (m 2 H /Q 2 ) remainders are provided by fitting the values 
to analytic functions with a number of free parameters. These final expressions are slightly 
approximate, but the error in F^(x, Q 2 ) is only ~ 1% even in the most extreme cases. 

There is one more problem in defining the VFNS. The ordering for F^(x,Q 2 ) is different 
for the rif and rif + 1 regions. This can be illustrated by the following table which shows the 
order by order expressions both below and above the transition point. 



rif — flavour rif + 1 — flavour 

LO ^C FF / g ®g^ C v 2 F / H ®{h + h) 

\ 2 



NLO (2a) (<w ® g»'+CZ% ® (C 2 ^ ® (/i + fc) + C7j# (g> 

The issue is that the series expansion begins at zeroth order above the transition point, where 
there is a heavy flavour distribution, but at 0(as) below the transition point. Hence, what 
is meant by LO, NLO etc. is different by one power of «s as one changes the number of 
active quark flavours. Therefore, making the transition directly from a given fixed order to the 
same relative order when going from rif to rif + 1 flavours leads to a different order in as and 
discontinuities which may be rather significant- for example, LO is nonzero as one approaches 
the transition point from below, but zero when approaching it from above. One must make 
some decision on how to deal with this problem. 

Up to now ACOT have used the same order of 0:5 above and below the transition point, 
e.g. at NLO 



^Cl F / 9 ® g»f - Cl F d l ®(h + h) + 2( C Ihh ®(h + h) + C FF / g ® g^ 1 ). (21) 



The structure function is then automatically continuous. However, there is effectively LO 
evolution below the transition point - C F Jj g contains only information on Pjj g , not on P 1 - and 
NLO evolution above it. Hence the slope dF^(x,Q 2 )/d \nQ 2 is discontinuous. 

The Thorne-Roberts scheme used the same relative order above and below the transition 
point, but added a uniquely determined (^-independent term above the transition point to 
maintain continuity of the structure function. For example, at LO 

^^C FF / g {Q 2 /rn 2 H )®g^{Q 2 ) - ^^^£(1)®^ 

(22) 

i.e. this prescription freezes the higher order 0:5 term when going upwards through Q 2 = m 2 H . 
This difference in choice is extremely important at low Q 2 (if using /x 2 = Q 2 ), as is illustrated in 
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Fj TR style and ACOT style 




-1 -4 -3 -2 -I -"5 -4 -3 -2 

10 10 10 10 10 1 10 10 10 x 10 10 



Figure 2: Comparison between the ACOT choice of ordering and the Thorne-Roberts choice at 
NLO. 
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Fig. 2 which compares the two choices at NLO. The 0(a 2 s ) part is dominant at low x for Q 2 < m 2 c 
because the 0(a 2 s ) coefficient functions diverge at small x whereas the 0(as) coefficient function 
is finite in this limit. Indeed, the "frozen" part is very significant for m 2 < Q 2 < 12GeV 2 . 
Its inclusion clearly improves the match to the data [12, 13]. It is also clear that switching 
from the standard n^-flavour NLO to the standard rif + 1-flavour NLO would lead to a large 
discontinuity in F 2 H (x,Q 2 ). Hence we choose to continue using the Thorne-Roberts approach 
for the matching at the transition point. However, this is a place where there is a definite 
freedom of choice, and there are various possibilities available. 

With the type of choice made for the definition of the heavy flavour coefficient functions 
C^[^Ij{Q 2 /m 2 H ,z) in Eq.(17) there is no problem with ordering across the transition point for 
the longitudinal structure function. This is because both the FFNS and VFNS coefficient 
functions begin at 0(a s ) for both the gluon and heavy quarks. It is possible to choose a 
non-vanishing value for C^° H (Q 2 /m 2 H , z), and indeed some VFNS definitions do so, but this 
is contrary to the spirit of this approach, and will lead to extra complications. 

In order to define fully the VFNS at NNLO, this choice for ordering above and below the 
transition point means that we need the O{o i s ) heavy-flavour coefficient functions for Q 2 < m 2 H 
and that the contribution from these should be frozen for Q 2 > m 2 H . However, these coefficient 
functions are not yet known. Nevertheless, we do know the leading threshold logarithms [14], 
i.e. the leading contribution for W 2 not much above Am 2 H . This is given by 

C 2 , Hg (Q /m H , z) ~ z{1+r])Q 2 + Am 2 H P(v, Q /m H ), V = ~ 1, (23) 

i.e. rj — > at threshold and rj — > oo as W 2 — >• oo. p(r] : Q 2 /m 2 H ) is a function which models 
the contribution from the dominant threshold logarithms. This contribution occurs only in the 
gluon sector. 

We can also derive the leading ln(l/x) term from /^-dependent impact factors derived by 
Catani, Ciafaloni and Hautmann [15] . With a little work these can be shown to give 

C^t l ° WX (Q 2 /m 2 H , z) = 96^^ 2 (Q 2 /<), (24) 

where ^{Q 2 /m 2 H ) may be calculated and ^(l) ~ 4. We also know that in this small-x limit 
C2,Hq lOWX (Q 2 /m 2 H ,z) = 4/9C^f lowx (Q 2 /m 2 H ,z). By analogy with the known NNLO coeffi- 
cient functions and splitting functions it is reasonable to propose that this be modified to 

C^f l ° wx (Q 2 /m 2 H , z) = -/3(ln(l/z) - 4)(1 - z/x max ) 20 <k 2 {Q 2 /m 2 H ). (25) 

z 

(3 = (1 — Am 2 H z I \Q 2 {\ — z))) 1 / 2 is the velocity of a heavy quark in the centre-of-mass frame 
and its introduction ensures that this contribution — > smoothly at threshold. The leading 
ln(l/ z) is accompanied by ~ —4, i.e. al/z term of similar size to that in other known coefficient 
functions and splitting functions is introduced. Finally, the effect of the this entire small z term 
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Figure 4: Comparison of the NLO and NNLO predictions for F£(x, Q 2 ) compared with HI 
data. 
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is damped as z — > 1 by the large power of (1 — z/x max ). The power of 20 is chosen in order to 
make the contribution in Eq.(25) very suppressed until x < 0.1, which is in line with the values 
of x above which the small- a; divergent terms tend to be suppressed in the complete NNLO 
splitting functions [20]. The total approximate NNLO coefficient function is obtained by adding 
the contributions in Eq.(23) and Eq.(25) The amount of information is similar to (though a 
little weaker than) that used previously to derive approximate NNLO splitting functions [21], 
which turned out to be a very good approximation once the exact expressions became known. 

These expressions could also be used to provide an approximate NNLO FFNS definition. 
However, in this case they would have to be used over a far wider range of Q 2 , rather than the 
small range here. In particular the frozen NNLO contribution becomes a proportionally very 
small contribution to the total F^(x,Q 2 ) at high Q 2 , whereas there is no reason to believe 
the 0(a 3 s ) contribution in the FFNS is small since it contains terms of order \n 3 (Q 2 /m 2 H ). 
The NNLO FFNS would therefore be more genuinely approximate. For F^(x,Q 2 ) we can 
approximate the 0(a 3 s ) coefficient functions at Q 2 < m 2 H using the same approach. In this case 
there is no large threshold contribution and the small-x contribution is estimated to be 

C ™¥™°(C?/m 2 B , z) = -(3 3 (\n(l/z) - 4)(1 - z/x max ) 20 K L (Q 2 /m 2 H ), (26) 

z 

where ~ 0.16 and the (3 3 reflects the fact that the longitudinal heavy-flavour coefficients 

are much more suppressed near threshold. For the VFNS for F^(x, Q 2 ) one can then extrapolate 
smoothly to the exact massless 0(a 3 s ) coefficient functions used in the VFNS at high Q 2 . 
Explicitly, for Q 2 > m 2 H 
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i n f 



+ Ktt4^4)(^) 3 ?^ 3 ^"' +1 (« 2 )' ™ 

where the term oc C^ l g ® A 2 g H eg) g nf maintains continuity of the structure function across 
the transition point despite the discontinuity in the gluon distribution. Again the degree of 
modelling and approximation is far less than in a NNLO FFNS. 

Using these approximate 0(a 3 s ) coefficient functions one can produce full NNLO predictions 
for structure functions with discontinuous partons and coefficient functions but continuous 
F H (x, Q 2 ). The results are not very sensitive to the choices made in this approximation, as long 
as they are within a sensible range. Note also that the definition of the VFNS, relying only on 
Eq.(9), Eq.(14), Eq.(17) and the ordering across the transition point, may be straightforwardly 
generalised to any choice of factorization and renormalization scales. For the simplest choice of 
/i 2 = Q 2 the NNLO corrections are seen in Fig.3. They clearly improve the match to the lowest 
Q 2 data, where NLO is always too low. This large increase at low x is because the NNLO 
coefficient functions are more divergent than the NLO coefficient functions. The comparison 
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with the recent bottom quark production data [16] from HI is also shown in Fig. 4. 4 The 
agreement is good, but one can clearly see that the slope (dF^x, Q 2 ) / d\n.Q 2 ) is smaller at 
NNLO than at NLO, and the same is true for charm. This is because one is interpolating from 
a higher value at low Q 2 , due to the NNLO coefficient function contribution, to lower values at 
high Q 2 , where the quark content is dominant and has a relative reduction compared to NLO 
because the evolution has started from a negative value rather than zero. Hence, this tendency 
for a reduction in the slope of the heavy-flavour structure function at NNLO is a generic feature 
of NNLO. The detailed phenomenology of the global fit with the NNLO VFNS prescription 
will appear in the account our next global parton analysis [18]. 



3 Light-Flavour Sector 

At NNLO it is no longer possible to think of heavy-flavour effects as only associated with heavy 
flavour production, as has already been discussed in [19]. At this order we also get contributions 
due to heavy flavours away from the photon vertex. Examples of such contributions are shown 
below. 




In the light-flavour sector the VFNS is defined as before. There are matrix elements changing 
the light flavours as we go from nj flavours to nf + 1 flavours and contributions to light-flavour 
coefficient functions in the FFNS from the type of diagrams above. The exact expressions for 
these, as well as the asymptotic limits, can be found in [10]. Defining the consistent VFNS 
according to Eq. (9) then leads to a discontinuity in the coefficient functions across the transition 

4 This is an updated version of Fig. 8 in [16]. This previous figure used the MRST04 NNLO partons which 
applied a rather approximate NNLO treatment of heavy flavours. The updated figure is constructed using 
NNLO partons from a fit which applies the NNLO VFNS. This leads to a generally slightly increased prediction 
for Q 2 ). The gluons from this new fit are exhibited in Fig. 3 of [17], and are more negative than our 

previous, approximate NNLO partons at low Q 2 and very small x. 
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point, which again cancels that in the light quark distributions, leaving the total structure 
function continuous. In this sector there are peculiar complications due to the occurrence of 
(ln m (l — z)/(l — z)) + terms at the threshold. For example, from the type of diagrams on the 
right-hand side we get contributions to the structure function of the form 

al\n\Q*/ml) ^\- z Z) )® q {x^) = a| ln 3 (Q 2 /m 2 H )q(x,Q 2 ) 

W s \n 2 (Q 2 /m 2 H ) ( ln 7^_~ Z) ) ® Q 2 ) + 0(\n 2 (Q 2 /m 2 H )). (28) 

The maximum value of z in the convolution is 1/(1 + 4m 2 H /Q 2 ) so the divergence at z = 1 
is not reached. However, the convolution on the left-hand side does produce an additional 
power of \n(Q 2 /m 2 H ). The \n 3 (Q 2 /m 2 H ) contributions then cancel exactly with contributions 
from terms of the form a 2 \v?(Q 2 /m 2 H )5(l — z) in the coefficient functions from the type of 
diagrams on the left-hand side, leaving the remainder, including the "-(-"-distributions. This 
is an added complication compared to the case where the heavy quark appears at the photon 
vertex, requiring particular care in a numerical implementation of the VFNS. However, it is 
only really a technical problem rather than producing any fundamentally new features. In 
the light-flavour sector there is no problem with ordering. The light-flavour contribution to 
F 2 (x, Q 2 ) begins at zeroth order in the FFNS and VFNS. 

Furthermore, if one is being totally correct, the left-hand type diagram and the soft parts 
of the right-hand type diagram should contribute to the light-flavour structure function, and 
the hard part of the right-hand type diagram contributes to F^(x,Q 2 ) [19]. This can be 
implemented (it depends on a separation parameter, determining "hard" and "soft"), but each 
contribution is in practice tiny. At the moment we include all such contributions in the light 
flavours. This leads to a very small underestimate of the heavy-flavour structure functions. 

4 Charged- Current Structure Functions 

The VFNS works, in principle, in much the same way for charged currents. The zeroth-order 
coefficient function for single (anti) charm production from a (anti) strange quark is now 

Cl!f{Q 2 /ml z) = S(z - Q 2 /(Q 2 + m 2 c )), (29) 

i.e. the threshold is now for a single heavy quark production. The same is obviously true for 
(anti)charm production from a (anti)down quark. This is far simpler than the prescription in [8]. 
At higher orders the generalization is again trivial. Using the ACOT(x) reasoning one simply 
uses the light coefficient function with argument z replaced by zj(\ + m 2 / Q 2 ). However, there 
is a major problem in defining the full VFNS, even at NLO. The FFNS coefficient functions are 
calculated at 0(ots) [22] but are not yet calculated exactly at 0(a 2 s ) for the charged-current 
case. 5 The generalization from the neutral current case is not simple because the divergences 

5 They have been determined in the Q 2 ^> m 2 H limit in [23]. 



13 



in the final state are different, i.e. one more particle is massless and regularized by dimensional 
regularization, and one fewer is regularized by the mass. 

Hence, in defining the VFNS for the charged current case it is necessary to make an ap- 
proximation. In practice we have used the 0{a 2 s ) neutral current cross-sections but altered the 
threshold dependence in all expressions so that all dependence on Q 2 /(Q 2 + Am 2 H ) is replaced 
by dependence on Q 2 /(Q 2 + rn 2 H ). This guarantees that all terms respect the true kinematic 
threshold. The approximation occurs mainly at low Q 2 where the 0{a 2 s ) FFNS coefficient 
functions are most important in the VFNS. This is not peculiar to this definition of a VFNS, 
but will be present in any current approach. 6 There is an uncertainty at low Q 2 at 0(a 2 s ) that 
can only be removed by an explicit calculation. At higher Q 2 the VFNS tends to the zero-mass 
limit for all coefficient functions, so all expressions become exact. Hence the HERA charged 
current data are very insensitive indeed to the approximation. The charged current data at 
low Q 2 only exist down to i ~ 0.01, much higher than for the neutral current HERA data, 
and at low Q 2 and high-x heavy-flavour production is very small. Hence any errors in the 
approximation are not very important phenomenologically. 

5 Conclusions 

There are discontinuities in both the parton distributions and the coefficient functions at NNLO. 
This makes a variable-flavour number scheme more necessary than ever. The ZMVFNS is badly 
discontinuous at the transition point Q 2 = m 2 H) and the FFNS is only approximate at NNLO. A 
generalization of the ACOT(x) prescription leads to a physically sensible and simple VFNS, in 
principle defined to all orders. One must still be careful about matching when going across the 
transition point of Q 2 = m 2 H . If this matching is done properly it guarantees the continuity of 
the physical structure functions and maximises the smoothness of the function. We choose the 
Thorne-Roberts method of matching above and below the transition, i.e. choose the correct 
order for the region of n/ flavours and add an additional, uniquely defined constant for the 
region of rif + 1 flavours to guarantee continuity. This choice is significant and leads to a much 
better match to the low Q 2 data. We have devised an explicit, full NNLO VFNS for F 2 (x, Q 2 ) 
and Fl(x,Q 2 ), with a small amount of necessary modelling of NNLO fixed-flavour coefficient 
functions. The NNLO variable-flavour number scheme seems to improve the fit to lowest x and 
Q 2 data greatly and is not very sensitive to this modelling. It is essential to use such an NNLO 
VFNS in NNLO global analysis of data, and indeed the construction of the NNLO VFNS makes 
such a precise analysis possible [18]. 

6 This is except for the ZMVFNS, which will simply be completely wrong by terms of 0(ml/Q 2 ) and will 
have incorrect kinematic thresholds, whereas the approximation here is guaranteed to have the correct general 
form. 



14 



Acknowledgements 

I would like to thank Stefan Kretzer, Alan Martin, Fred Olness, James Stirling, Wu-Ki Tung 
and Chris White for various discussions on this subject. I would also like to thank Mandy 
Cooper-Sarkar, Claire Gwenlan and Paul Thompson for numerous discussions on the compar- 
ison to data, and the latter for providing me with Fig.4. I would also like to thank the Royal 
Society for the award of a University Research Fellowship. 

References 

[1] E. Laenen, S. Riemersma, J. Smith and W. L. van Neerven, Nucl. Phys. B392 (1993) 
162; 

B. W. Harris and J. Smith, Nucl. Phys. B452 (1995) 109. 
[2] J. C. Collins and W. K. Tung, Nucl. Phys. B278 (1986) 934. 
[3] M. Buza, et ai, Eur. Phys. J. CI (1998) 301. 

[4] R.S. Thorne, proceedings of 13th International Workshop on Deep Inelastic Scattering 
(DIS 05), Madison, Wisconsin, 27 Apr - 1 May 2005, p. 847, hep-ph/0506251. 

[5] J. C. Collins, Phys. Rev. D58 (1998) 094002. 

[6] M. Aivazis, , J. C. Collins, F. Olness and W. K. Tung, Phys. Rev. D50 (1994) 3102. 

[7] R. S. Thorne and R. G. Roberts, Phys. Lett, bf B421 (1998) 303; Phys. Rev. D57 (1998) 
6871. 

[8] R. S. Thorne and R. G. Roberts, Eur. Phys. J. C19 (2001) 339. 

[9] W. K. Tung et ai, J. Phys. G28 (2002) 983; S. Kretzer et ai, Phys. Rev. D69 (2004) 
114005. 

[10] M. Buza, et ai, Nucl. Phys. B472 (1996) 611. 

[11] S. Riemersma, J. Smith and W.L. van Neerven, Phys. Lett. B347 143, (1995). 

[12] ZEUS collaboration: S. Chekanov et ai, Phys. Rev. D69 (2004) 012004. 

[13] HI collaboration: C. Adloff et ai, Phys. Lett. B528 (2002) 199. 

[14] E. Laenen and S. Moch, Phys. Rev. bf D59 (1999) 034027. 

[15] S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. bf B366 (1991) 135. 

[16] HI collaboration: A. Aktas et ai, Eur. Phys. J. C45 (2006) 23. 



15 



[17] R. S. Thorne, the proceedings of "New Trends in HERA Physics 2005" Ringberg Castle, 
Tegernsee, October 2005, p. 359, hep-ph/0511351. 

[18] A. D. Martin, W. J. Stirling and R. S. Thorne, to appear. 

[19] A. Chuvakin, J. Smith and W. L. van Neerven, Phys. Rev. D61 (2000) 096004. 

[20] A. Vogt, S. Moch and J.A.M. Vermaseren, Nucl. Phys. B688 (2004) 101; Nucl. Phys. 
B691 (2004) 129. 

[21] W. L. van Neerven and A. Vogt, Nucl. Phys. B568 (2000) 263; Nucl. Phys. B588 (2000) 
345; Phys. Lett. B490 (2000) 111. 

[22] M. Gliick, S. Kretzer and E. Reya, Phys. Lett. B B380 (1996) 171. 

[23] M. Buza and W. L. van Neerven, Nucl. Phys. B500 (1997) 301. 



16 



